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TWO ANTENNA, TWO PASS INTERFEROMETRIC SYNTHETIC APERTURE 

RADAR 

This invention was developed under Contract DE-AC04-94AL85000 between Sandia 
5 Corporation and the U.S. Department of Energy. The U.S. Government has certain rights in 

this invention. 

FIELD OF THE INVENTION 

The invention relates generally to data collection for topographic modeling and, more 
10 particularly, to two antenna, two pass interferometric synthetic aperture radar ("IFSAR"). 

BACKGROUND OF THE INVENTION 

Interferometric synthetic aperture radar ("IFSAR") is a radar system used to obtain 
target height information and form multi-dimensional maps of imaged areas. IFSAR utilizes 
15 at least two (2) SAR images of the same scene, formed at slightly different elevation angles 

relative to each other, to extract information about target heights. Such images can be 
coherently combined to ascertain the topography of the imaged area and produce three- 
dimensional maps of the imaged area. 
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Currently, IFSAR utilizes two (2) principal operational modes. The first mode is 
two-antenna, one-pass IFSAR in which a single aircraft with two (2) antennas, displaced in a 
direction normal to the flight path of the aircraft, flies by the scene of interest once while 
collecting data. The second mode is one-antenna, two-pass IFSAR in which a single aircraft 
with a single anteima flies by the scene of interest twice, along slightly offset flight paths, 
while collecting data. 

Each of these IFSAR operation modes possesses strengths and weaknesses. The two- 
antenna, one-pass mode allows for precise knowledge of baseline length, but often limits that 
baseline length. Additionally, residual system noise manifests itself as uncertainty in the 
topographic height estimates with lesser baseline lengths resulting in increased height-noise 
and limiting the effective topographic resolution. Because the one-antenna, two-pass mode 
allows a greater separation of the anteimas, it allows a larger baseline, and thereby has a 
diminished sensitivity to system noise. Although it produces less height-noise, it suffers 
from imprecise knowledge of the baseline length (antenna separation) between passes, and 
consequently does not provide accurate target height scaling. Conventional navigation 
instruments are inadequate in measuring aircraft flight paths with suitable precision. 
Although a two-antenna, two-pass IFSAR has the potential to provide both the scaling 
accuracy of the two-antenna, one-pass mode and the height-noise performance of the one- 
antenna, two-pass mode, it currently does not have precise knowledge of the baseline length 
between passes. 
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It is well-known by workers in the art that the uncertainty in the target height 
decreases as the antenna baselines increases. By decreasing the uncertainty in the target's 
height, the accuracy of the resulting digital elevation model improves. One alternative for 
increasing the antenna baseline length entails equipping a single aircraft with additional 
5 antennas with suitable offsets to generate greater baselines.. This is expensive and 

complicates the engineering. Additionally, flexing of the aircraft structure during flight may 
also be severe enough to render unacceptable uncertainty in the larger baselines without 
somewhat expensive baseline measurement schemes or devices. 

It is therefore desirable to provide a solution that avoids the weaknesses of the 
10 aforementioned conventional IFSAR modes. 

Exemplary embodiments of the invention permit data driven alignment of multiple 
independent passes, thereby providing the scaling accuracy of the two-antenna, one-pass 
mode and the height-noise performance of the one-antenna, two-pass mode. The antenna 
baseline between multiple flight passes can be accurately estimated from the data itself, 
15 thereby reducing both height-noise and scaling error and allowing for more accurate 

information about target ground position locations and heights. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The above and further advantages of the invention may be better understood by 
referring to the following description in conjunction with the accompanying drawings in 
which corresponding numerals in the different figures refer to the corresponding parts, in 
which: 

FIGURE 1 diagrammatically illustrates exemplary flight passes for a two-antenna, 
two-pass IFSAR system in accordance with exemplary embodiments of the present 
invention; 

FIGURE 1 A displays a table of some exemplary parameters for a two-antenna, two- 
pass IFSAR system in accordance with the present invention; 

FIGURE 2 diagrammatically illustrates geometry for a two-antenna, one-pass IFSAR 
system in accordance with the knpwn art; 

FIGURE 3 diagrammatically illustrates geometry for a two-antenna, one-pass IFSAR 
system with a negative antenna angle in accordance with the known art; 

FIGURE 4 diagrammatically illustrates geometry for a two-antenna, two-pass IFSAR 
system in accordance with exemplary embodiments of the present invention; 

FIGURE 5 diagrammatically illustrates two antenna pairs and pertinent geometry for 
an IFSAR system in accordance with exemplary embodiments of the present invention; 

FIGURE 6 diagrammatically illustrates a single antenna pair and pertinent geometry 
for an IFSAR system in accordance with the known art; 

-4- 



3386467vl 



Patent Application 
Docket No. SD-7196 

FIGURE 7 graphically illustrates RMS values for baseline uncertainty in accordance 
with exemplary embodiments of the present invention; 

FIGURE 8 graphically illustrates RMS uncertainty values for a target's ground range 
position in accordance with exemplary embodiments of the present invention; 
5 FIGURE 9 graphically illustrates RMS uncertainty values for a target's height in 

accordance with exemplary embodiments of the present invention; and 

FIGURE 1 0 diagrammatically illustrates an exemplary embodiment of an apparatus 
for implementing a two-antenna, two-pass IFS AR in accordance with the present invention. 
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DETAILED DESCRIPTION 

While the making and using of various embodiments of the present invention are 
discussed herein in terms of two-antenna, two-pass interferometric synthetic aperture radar 
("IFSAR"), it should be appreciated that the present invention provides many inventive 
concepts that can be embodied in a wide variety of contexts. The specific embodiments 
discussed herein are merely illustrative of specific ways to make and use the invention, and 
are not meant to limit the scope of the invention. 

The present invention provides a solution that enables data driven alignment of 
multiple independent passes, thereby providing the scaling accuracy of the two-antenna, one- 
pass IFSAR mode and the height-noise performance of the one-antenna, two-pass IFSAR 
mode. The two-antenna, two-pass IFSAR mode presented herein can accurately estimate the 
antenna baseline between two (2) flight passes from the data itself and reduce height-noise 
and scaling error, allowing for more accurate information about target ground position 
locations and heights. The two-antenna, two-pass IFSAR mode presented herein can use the 
data collected from each individual flight pass to estimate the antenna baseline length 
between the passes. True three-dimensional radar imaging fi-om stand-off aircraft and 
satellite platforms can be obtained via multi-pass IFSAR which can allow for data collectiom 
from more than two (2) flight passes. 

FIGURE 1 diagrammatically illustrates exemplary flight passes 1 and 2 for a two- 
antenna (20a and 20b), two-pass (1 and 2) IFSAR system in accordance with exemplary 
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embodiments of the present invention. Aircraft 10, carrying antennas 20a and 20b 
(alternatively denoted as antenna pair 20), can perform flight passes 1 and 2 along flight 
trajectory 30 over scene of interest 40. Any type platform can be used to carry the radar 
system. Each flight pass can be performed at a different range, such as Rl and R2, from 
5 scene of interest 40. In some embodiments, each flight pass, such as flights passes 1 and 2, 

can possess a different elevation angle to scene of interest 40. Flight passes 1 and 2 are 
separated by a distance of B meters, therefore B represents the antenna baseline length 
between antenna pair 20 during flight pass 1 and antenna pair 20 during flight pass 2. In 
some embodiments, B may be on the order of tens of meters. In some embodiments, flight 

10 passes 1 and 2 are as close to parallel as possible (similar to the one-antenna, two-pass 

mode). However, since IFSAR has been shown to be possible with non-parallel flight paths, 
the data itself can also be used to calibrate the relative flight paths. Some exemplary 
parameters for a two-antenna, two-pass IFSAR system are given in the table of FIGURE 1 A. 
Despite modem navigation equipment, the antenna baseline length between two (2) flight 

1 5 passes is generally unknown at an accuracy and precision required to produce accurate target 

ground range position and height infonriation for the imaged scene. While IFSAR 
topography measiu-ements may be calibrated using ground-truth measurements, opportunities 
to collect such ground-truth measurements are not always available. Consequently, the 
antenna baseline length between flight passes needs to be accurately estimated by altemate 

20 means to yield improved ground range position and height accuracies. 
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FIGURE 2 diagrammatically illustrates geometry for a two-antenna, one-pass IFS AR 
system in accordance with the known art. In FIGURE 2, the center of the scene is denoted as 
the scene's central reference point ("CRP"). In practice, distances (such as ri and r2) 
between the target (such as target 45) and antennas (such as antennas 20a and 20b) are very 
large relative to the smaller antenna baseline lengths (such as b) and scene sizes (such as 
scene of interest 40, FIGURE 1). Under these conditions, the equations for estimating the 
target's ground range position Sy and height Sz when antenna angle 9l is positive are: 

on1/2 
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where \|/c represents the grazing angle of the radar system, which is conventionally 
available from the navigational system on the aircraft. Ranges ri and rc are also 
conventionally available from the navigational system, while distance s between the CRP 
and target 49 can be calculated from data conventionally available from the navigational 
system as discussed below with reference to FIGURE 4. 

Figure 3 displays the geometry for a two-antenna, one-pass IFSAR system when 
antenna angle 9l is negative. Range rz and angle y are conventionally available from the 
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aircraft's navigational system and distance s between the CRP and target 49 can be 
calculated from data conventionally available from the navigational system. Using the law 
of cosines, ^ is related to target angle 9t by (n - 6t) in radians or (180^ - 9t) in degrees. 
Target angle Qj can be calculated as shown below: 



(3) 



Using these equations, the target's ground range position Sy and height estimation 
equations when antenna angle 9l is negative are: 
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Equations 1 and 2 use range ri with range rc, while equations 4 and 5 use range r: 
with range rc due to the setup and geometry of the problem. The long range condition 
between the platform and the scene, the small scene size, and the small baseline length 
condition between the two (2) individual antennas are all valid, making range ri 
approximately equal to range r2, both of which are conventionally available from the 
aircraft's navigational system. Consequently, equations 1 and 2 can be interchanged with 
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equations 4 and 5. The explicit derivations for these equations are well known to workers in 
the art. 

FIGURE 4 diagrammatically illustrates geometry for a two-antenna, two-pass EFS AR 
system in accordance with exemplary embodiments of the present invention. Antennas 11 
5 and 12 comprise an antenna pair (alternatively designated as 100) for a first flight pass and 

are separated by a distance of bi meters. Antennas 21 and 22 comprise an antenna pair 
(alternatively designated as 200) for a second flight pass and are separated by a distance of 
b2 meters. The distances separating the antennas in each pair (i.e., bi and b2) may differ and 
may be on the order of submeters (e.g., 0.3 meters). 

1 0 Since it is conventionally known that the root mean square ("RMS") error associated 

with a target's ground range position and height estimates are inversely proportional to the 
antenna baseline length, it follows that a larger baseline length will yield smaller errors for 
these estimates. Consequently, by utilizing the antenna baseline length B between two (2) 
passes rather than the antenna baseline length (i.e., bi or bi) between two (2) antennas of a 

15 single pass (see FIGURE 4) in conjunction with information from the flight passes, target 

ground range position Sy and height scaling accuracies can be improved to the equivalent 
of the two-antenna, one-pass IFSAR system while retaining the height-noise performance 
associated with the one-antenna, two-pass IFSAR system. For a two-antenna, two-pass 
mode, midpoint Midi between antennas 11 and 12 and midpoint Midi between antennas 21 

20 and 22 can be utilized as shown in FIGURE 4. The equations below can calculate an 
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estimate of the antenna baseline length B between midpoint Midi of antenna pair 100 and 
midpoint Midi of antenna pair 200. The distance from the scene's CRP to target 45 itself is: 

s = ^/sf+sT" . (6) 
Using conventional trigonometric equations and identities with the recorded phase 
(v|/ci and \|/c2) measurements, the ground range position and height of Midi are: 

and 

^z^ =rc^ sin(v^ci) • (8) 
Similarly, the ground range position and height of Midi are: 

=rc2COs(v^c2) (9) 

and 

=rc2-sin(<^c2) • (10) 
The range (rd and Vd) and phase (vj/d and \\fcz) measurements are conventionally 

available from the aircraft's navigational system. 

Using these two coordinate pairs, the estimation equation for antenna baseline length 

B between the two (2) flight passes is: 

B = V(Sz.-^)'+(Sy.-Sy,)'- (11) 

This result can be used to re-estimate the target's ground range position Sy and height 

Sz for a two-antenna, two-pass IFSAR system, resulting in refined estimates of the target's 

-11- 

3386467vl 



Patent Application 
Docket No. SD-7196 

ground range position Sy and height Sz. When using multiple targets to calculate the antenna 
baseline estimate, equations 7-11 can be calculated for each of the targets and the resulting 

A 

baseline estimates averaged to obtain B. FIGURE 5 is similar to FIGURE 4 and 
diagrammatically illustrates two antenna pairs (antennas 11 and 12, alternatively designated 
as antenna pair 100, and antennas 21 and 22, alternatively designated as antenna pair 200) 
and pertinent geometry for an IFS AR system in accordance with exemplary embodiments of 

A 

the present invention. Using antenna baseline length B between the two (2) flight passes 
(equation 11), the equation for antenna angle 0l becomes: 



^L=sin"* 



B ^ A((» + 2;r'm) ?i}{^^'lK'mf 
4;rB 32;r'Br2 



(12) 



This equation can account for phase ambiguities where the variable 'm' represents 
the integer ambiguity index. As conventionally known, there are ambiguities in the phase 
due to the fact that the phase is equivalent every 2;r radians. The possible values of 'm' are 
bounded by the navigational system and can be generated as conventionally known to 
workers in the art. In some embodiments, equation 12 can be calculated once for each *m' 
value. Also, A is the system wavelength (e.g., 0.02 meters as in FIGURE lA), and ^is the 
phase difference conventionally calculated as the difference in the range times the constant 
2;r/yl(i.e., ^)=(2;r//l)*(r2-ri)). 
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If equation 12 yields a positive antenna angle 0l, and noting that ranges and Fc are 
conventionally available from the aircraft's navigational system while distance s from the 
scene's CRP to target 45 can be calculated from data conventionally available from the 
navigational system, the following refined estimation equations apply: 
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If antenna angle Gt is negative in equation 12, the appropriate refined estimation 
equations become: 

1I/2 
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Equations 13 and 14 are similar to equations 1 and 2, but use range rz rather than 
range ri, while equations IS and 16 are the same as equations 4 and 5 (all use range r2). As 
discussed above with reference to equations 1, 2, 4, and 5, the long range condition between 
the platform and the scene, the small scene size, and the small baseline length condition 
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between the two (2) individual antennas are all valid, making range ri approximately equal 
to range ri. Also in equations 13-16, can be either td or rc2 with its corresponding V|;c 
value (v|/ci or \|/g2, respectively). (It is and \|/c2 in FIGURE 10 described below.) 

FIGURE 10 diagrammatically illustrates exemplary embodiments of an apparatus for 
implementing a two-antenna, two-pass IFS AR in accordance with the present invention. The 
aircraft's navigational system 1005 can supply ranges ri, r2, rd and r^z^ phase measurements 
\|/ci and V|;c2, target ground range Sy, target height Sz, and system wavelength k. With rd, rc2, 
\|/d, and \|;c2 as inputs, ground range estimator 1010 implements equations 7 and 9, 
outputting ground range estimates Syi and Sy2, and height estimator 1015 implements 
equations 8 and 10, outputting height estimates Szi and Sz2. Baseline estimator 1020 receives 
the outputs of ground range estimator 1010 and height estimator 1015 and, implementing 

A 

equation 1 1 , produces antenna baseline estimate B which is input to antenna angle estimator 
1030. With \|/d and \|/c2 as inputs, *m' generator 1040 generates a range of 'm' values which 
are input to antenna angle estimator 1030. Phase difference estimator 1025 receives ranges 
ri and ri and system wavelength X from navigational system 1005. Phase difference 
estimator 1025 implements the equation ((Itt/X) * (r2 ~ ri)) to produce ^ which is input to 

A 

antenna angle estimator 1030. Using B , m, and ^ along with X and Vi from navigational 
system 1005, antenna angle estimator 1030 implements equation 12 and produces an 
estimate of antenna angle Gl- Selector 1035 receives from antenna angle estimator 1030 a 
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positive/negative indication for 9l. Based on the sign of 9l, selector 1035 will select either 
ground range and height estimator 1045 or ground range and height estimator 1050. Ground 
range and height estimator 1045 implements equations 15 and 16 for a negative 9l and 
ground range and height estimator 1050 implements equations 13 and 14 for a positive 6l. 
5 Reference distance estimator 1055 receives Sy and Sz from navigational system 1005, 

implements equation 6, and produces distance s (the distance from the scene's central 
reference point to the target) which is input to both ground range and height estimators 1045 
and 1050. In the exemplary embodiments of FIGURE 10, both ground range and height 
estimators 1045 and 1050 receive ranges ri, rc2, and phase measurement vj/c2 from 
10 navigational system 1005 and distance s from reference distance estimator 1055 for use in 
producing refined ground range and height estimations Sy and s^ . 

The uncertainties in a target's ground range position and height estimates for a two- 
antenna, one-pass IFSAR system in the presence of noise can be derived with reference to 
FIGURE 6 wherein antennas 20a and 20b (i.e., a single antenna pair) are located on the 
1 5 order of kilometers away from a scene of interest that includes central reference point CRP. 

These uncertainties are known to workers in the art. In FIGURE 6, the scene size is small 
compared to range fc, and height Sz from the ground to target 45 on a first order 
approximation is equal to: 

re-^L-cos(i^c)- (17) 
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Here, grazing angle \|/c is constant and antenna angle 0l varies. Under small scene 
size and long range constraints, point Pi is approximately the same distance from the radar 
platform as the CRP on a first order basis. Therefore, the distance between CRP and Pi is 
very small in comparison to range Tc. Under small antenna angle conditions and using first 
order partial derivative concepts as known in the art, the uncertainties in target 45 ground 
range position Sy and height Sz for a two-antenna, one-pass IFSAR system reduce to the 
following well known approximations: 

®' 16bVSNR 

and 

16bVSNR 

respectively, where, as known by workers in the art, the phase variance is related to 

the signal to noise ratio ("SNR") as: 

2 1 

o-i = ■ . (20) 

^ SNR ^ ^ 

The uncertainties in a target's ground range position Sy and height estimates for 

exemplary embodiments of a two-antenna, two-pass IFSAR system in the presence of noise 

can be derived with reference to FIGURE 4. Using the antenna baseline length estimation 

equation (Equation 11), the imcertainty in antenna baseline length B can be calculated in 

terms of the imcertainty in the difference between the anterma pairs' midpoint (Midi ^d 
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Mid2) height coordinates (i.e., the difference Szi - S22) and the uncertainty in the difference 
between the ground range coordinates for these two midpoints (i.e., the difference Syi - Sy2). 
The overall antenna baseline uncertainty with respect to the uncertainties in both of these two 
(2) differences is: 



(21) 



After computing the two (2) individual components in the antenna baseline length 
uncertainty, the uncertainty in grazing angle \|/ci can be utilized to obtain the final baseline 
uncertainty equation. Using the geometry in Figure 4, this overall uncertainty becomes: 

1 



s^-M-SNR 



•{(r,,-r,^)^.cos^(V.j)- 



=2- 



16;r'-b, 



4s= 



(22) 



(-s^+r;;-r,^)^.r,^a^ 



This equation assumes the use of M independent targets to calculate the antenna 
baseline length estimate. That is, in some embodiments, baseline estimator 1020 of FIGURE 
10 can implement equations 7-11 for each of M targets and then average the M, results to 

A 

obtain its output baseline length estimate B . The baseline uncertainty in terms of the RMS 
from a single target's measurements is shown in Figure 7. The RMS for the antenna baseline 
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length uncertainty is plotted as a function of SNR and the distance from the target to the 
scene's CRP. 

Two observations can be made from the baseline standard deviation presented in 
Figure 7. First, as the SNR increases, the standard deviation of the baseline length estimate 
decreases because as the signal power increases relative to the noise power, the overall SNR 
increases. As signal noise decreases, the uncertainties in estimating the antenna baseline, the 
target location, and the target height also decrease. The second observation is that as the 
target moves farther away from the scene's CRP (i.e., the distance from the scene's CRP 
increases), the uncertainty in the baseline length estimate decreases because the baseline 
length imcertainty is inversely proportional to the target's distance from the scene's CRP. 
As the distance increases, the baseline length uncertainty will converge to zero since the first 
term in the brackets in Equation 22 will dominate the other two terms. From these two 
observations, if can be inferred that the exemplary conditions for accurately estimating the 
antenna baseline length between the two antenna pairs is to utilize a high SNR with targets 
located as far away as possible from the scene's CRP within the scene of interest, hi other 
words, large target swaths are better for estimating baseline lengths than narrow target 
swaths. 

In a two-antenna, two-pass IFSAR system, the refined estimates of target grotmd 
range position Sy and height Sz estimates also possess uncertainties. These uncertainties, 
respectively, are: 
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(23) 



and 



cr^^=(rcC0S(<^c))2c7^ 



(24) 



Using the antenna baseline length B between the two antenna pairs (such as 100 and 
200), the uncertainty in the antenna angle Gl becomes: 



1 



SNR I 4nB 16;r'Br. 



2 J 



1 



2 Jt2 > 



327r^ L 



.(25) 



2 J 



It is evident from the above three equations that as the antenna baseline length 
estimation uncertainty increases, the uncertainties in both the target ground range position Sy 
and height Sz also increase. Utilizing the uncertainty in antenna angle 6l in Equation 25 with 
the uncertainty in the antenna baseline length B in Equation 22, the uncertainty in estimating 
the target's refined ground range position Sy is: 
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SNR 
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(26) 

where M targets are used to estimate the antenna baseline length B between two 
antenna pairs (such as 100 and 200). On a first order approximation, the uncertainty in the 
target groimd range position Sy with respect to the antenna baseline length B between two 
antenna pairs (such as 100 and 200) becomes the well known expression: 



^2 ^ A^Tc^-sin^(t^e) 



(27) 



Figure 8 displays the RMS values for the target's refined ground range position Sy 
error uncertainty. In this figure, as the SNR increases, the RMS values in the uncertainty 
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decrease. This is also true for the target's distance from the scene's CRP. Targets located 
farther away from the scene's CRP within the scene of interest possess smaller RMS values 
in the ground range position Sy uncertainties than targets located closer to the scene's CRP. 
As seen from the first order approximation, the ground range position Sy uncertainty will 
converge to the value in Equation 27. 

Using M targets to estimate the antenna baseline length B between two antenna pairs 
(such as 100 and 200), the uncertainty in the target's refined height Sz estimate is: 

v2 
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(28) 
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On a first order approximation, the uncertainty in the target height with respect to 
the antenna baseline length B between two antenna pairs (such as 100 and 200) becomes the 
well known expression: 

16;r2.SNR.B2 • (29) 

It can therefore be inferred that in order to decrease the uncertainty in the target 
height Sz, bigger antenna baselines B are necessary. Similar to the groimd range position Sy 
uncertainty, the target height Sz uncertainty converges to the first order approximation shown 
in Equation 29. 

Figure 9 displays the RMS uncertainty values for estimating the target's refined 
height Sz in terms of the SNR and the target's distance from the scene's CRP. Similar to the 
observations from the target's refined ground range position Sy estimation procedure, the 
RMS uncertainty values for the refined height Sz estimates decrease as the SNR is increased 
and as the target moves farther away from the scene's CRP within the scene of interest. 

It will be evident to workers in the art that the exemplary embodiments described 
above can be readily implemented by suitable modifications in software, hardware or a 
combination of software and hardware in conventional topographic modeling systems, such 
as IFSAR systems. 

Although exemplary embodiments of the present invention have been described in 
detail, it will be understood by workers in the art that various modifications can be made 
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therein without departing from the spirit and scope of the invention as set forth in the 
appended claims. 
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